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Abstract 

The closed time-path (CTP) formalism is a powerful Green's function formu- 
lation to describe nonequilibrium phenomena in field theory and it leads to a 
complete nonequilibrium quantum kinetic theory. In this paper we make use 
of the CTP formalism to write down a set of quantum Boltzmann equations 
describing the local number density asymmetries of the particles involved in 
supersymmetric electroweak baryogenesis. These diffusion equations auto- 
matically and self-consistently incorporate the CP-violating sources which 
fuel baryogenesis when transport properties allow the CP-violating charges 
to diffuse in front of the bubble wall separating the broken from the unbroken 
phase at the electroweak phase transition. This is a significant improvement 
with respect to recent approaches where the CP-violating sources are in- 
serted by hand into the diffusion equations. Furthermore, the CP-violating 
sources and the particle number changing interactions manifest "memory" 
effects which are typical of the quantum transport theory and are not present 
in the classical approach. The slowdown of the relaxation processes may 
keep the system out of equilibrium for longer times and therefore enhance 
the final baryon asymmetry. We also stress that the classical approximation 
is not adequate to describe the quantum interference nature of CP-violation 
and that a quantum approach should be adopted to compute the sources 
since they are most easily built up by the transmission of low momentum 
particles. 
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1. Introduction and summary 

Because of the presence of unsuppressed baryon number violating processes at high 
temperatures, the Standard Model (SM) of weak interactions fulfills all the requirements 
for a successful generation of the baryon number at the electroweak scale |]J. The 
baryon number violating processes also impose severe constraints on models where the 
baryon asymmetry is created at energy scales much higher than the electroweak scale . 
Unfortunately, the electroweak phase transition is too weak in the SM ||. This means 
that the baryon asymmetry generated during the transition would be subsequently 
erased by unsuppressed sphaleron transitions in the broken phase. The most promising 
and well-motivated framework for electroweak baryogenesis beyond the SM seems to 
be supersymmetry (SUSY). Electroweak baryogenesis in the framework of the Minimal 
Supersymmetric Standard Model (MSSM) has attracted much attention in the past 
years with particular emphasis on the strength of the phase transition ||, [|, |(| and the 
mechanism of baryon number generation J7l |L EL TO. 



Recent analytical [IT], [L2| and lattice computations [Lj| have revealed that the phase 
transition can be sufficiently strongly first order if the ratio of the vacuum expectation 
values of the two neutral Higgses tan (3 is smaller than ~ 4. Moreover, taking into 
account all the experimental bounds as well as those coming from the requirement of 
avoiding dangerous color breaking minima, the lightest Higgs boson should be lighter 
than about 105 GeV, while the right-handed stop mass might be close to the present 



experimental bound and should be smaller than, or of order of, the top quark mass |T2 

Moreover, the MSSM contains additional sources of CP-violation besides the CKM 
matrix phase. These new phases are essential for the generation of the baryon number 
since large CP-violating sources may be locally induced by the passage of the bub- 
ble wall separating the broken from the unbroken phase during the electroweak phase 
transition. Baryogenesis is fueled when transport properties allow the CP-violating 
charges to efficiently diffuse in front of the advancing bubble wall where anomalous 
electroweak baryon violating processes are not suppressed. The new phases appear in 
the soft supersymmetry breaking parameters associated to the stop mixing angle and 
to the gaugino and neutralino mass matrices; large values of the stop mixing angle 
are, however, strongly restricted in order to preserve a sufficiently strong first order 
electroweak phase transition. Therefore, an acceptable baryon asymmetry from the 
stop sector may only be generated through a delicate balance between the values of 
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the different soft supersymmetry breaking parameters contributing to the stop mixing 
parameter, and their associated CP- violating phases ||. As a result, the contribution 
to the final baryon asymmetry from the stop sector turns out to be negligible. On 
the other hand, charginos and neutralinos may be responsible for the observed baryon 
asymmetry if the phase of the parameter /i is large enough |Q, ID . Yet, this is true 
within the MSSM. If the strength of the electroweak phase transition is enhanced by 
the presence of some new degrees of freedom beyond the ones contained in the MSSM, 
e.g. some extra standard model gauge singlets, light stops (predominantly the right- 
handed ones) and charginos/neutralinos are expected to give quantitatively the same 
contribution to the final baryon asymmetry. 

The baryon asymmetry has been usually computed using the following separate steps 

I, HI 

1 ) Look for those charges which are approximately conserved in the symmetric phase, 
so that they can efficiently diffuse in front of the bubble where baryon number violation is 
fast, and non-orthogonal to baryon number, so that the generation of a non-zero baryon 
charge is energetically favoured. Charges with these characteristics in the MSSM are the 
axial stop charge and the Higgsino charge, which may be produced from the interactions 
of squarks and charginos and/or neutralinos with the bubble wall, provided a source of 
CP-violation is present in these sectors. 

2) Compute the CP- violating currents of the plasma locally induced by the passage 
of the bubble wall. The methods present in the literature properly incorporate the 
decoherence effects which may have a crucial impact on the generation of the CP- 
violating observable. 

3) Write and solve a set of coupled differential diffusion equations for the local 
particle densities, including the CP- violating source terms derived from the computation 
of the current at step 2) and the particle number changing reactions. The solution to 
these equations gives a net baryon number which is produced in the symmetric phase 
and then transmitted into the interior of the bubbles of broken phase, where it is not 
wiped out if the first transition is strong enough. 

It is important to notice that the CP- violating sources are inserted into the diffusion 
equations by hand only after the CP- violating currents have been defined and computed. 
This procedure is certainly appropriate to describe the damping effects on the CP- 
violating observables originated by the plasma interactions, but does not incorporate 
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any relaxation time scale arising when diffusion and particle changing interactions are 
included (even though this approximation might be good if the diffusion time scales are 
larger than the damping time scales) and is theoretically not consistent. Furthermore, 
since a certain degree of arbitrariness is present in the way the CP- violating sources 
may be defined, different CP-violating sources have been adopted for the stop and the 
Higgsino sectors in the literature J7|, ||. This is certainly not an academic question since 
different sources may lead to different numerical results for the final baryon asymmetry, 
especially if the sources are expressed in terms of a different number of derivatives of 
the Higgs bubble wall profile and, therefore, in terms of different powers of the bubble 
wall velocity v w and bubble wall width L w . 

It is indisputable that one might be able to rigously derive a set of transport (diffu- 
sion) equations already incorporating the CP-violating sources in a self-consistent way 
only by means of a more complete treatment of the problem. The goal of this paper 
is to show that nonequilibrium Quantum Field Theory provides us with the necessary 
tools to write down a set of quantum Boltzmann equations (QBE's) describing the local 
particle densities and automatically incorporating the CP-violating sources. The ordi- 
nary quantum field theory at finite temperature is not useful to study the dynamics of 
particle densities. This is because we need their temporal evolution with definite ini- 
tial conditions and not simply the transition amplitude of particle reactions with fixed 
initial and final conditions. The most appropriate extension of the field theory to deal 
with these issues it to generalize the time contour of integration to a closed time-path 
(CTP). The CTP formalism is a powerful Green's function formulation for describing 
nonequilibrium phenomena in field theory, it leads to a complete nonequilibrium quan- 
tum kinetic theory approach and it will guide us towards the rigorous computation of 
the CP- violating sources for the stop and the Higgsino numbers. This will also eliminate 
the level of arbitrariness the previous treatments are suffering from. 

There exist other good reasons why one should call for the nonequilibrium quan- 
tum kinetic theory. The fact that CP-violating sources are most easily built up by the 
transmission of low momentum particles over a distance [7], |8|] is an indication that 
particles with masses smaller than or of the order of the temperature T are relevant in 
the process of quantum interference leading to CP-violating sources in the bubble wall. 
Basically, the sources are dominated by particles with long wavelengths in direction 
perpendicular to the wall. This tells us that the classical approximation is not adequate 
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to describe the quantum interference nature of CP- violation and a quantum approach 
must be adopted to compute the sources. For low momentum particles, the validity of 
the classical Boltzmann equation starts to break down and the ultimate answer can be 
provided only by a complete nonequilibrium quantum field theory approach. Kinetic 
theory and classical Boltzmann equations have been used to describe the dynamics 
of particles treated as classical with a defined position, energy and momentum. This 
requires that, in particular, the mean free path must be large compared to the Comp- 
ton wavelength of the underlying particle in order for the classical picture to be valid, 
which is not guaranteed for particles with a small momentum perpendicular to the wall. 
Distribution functions obeying the quantum Boltzmann equations are the only correct 
functions to describe particles in an interacting, many-particle environment. Further- 
more, we will show that the CP-violating sources and the particle number changing 
interactions built up from the CTP formalism are characterized by "memory" effects 
which are typical of the quantum transport theory |15, [TjJ. In the classical kinetic theory 



the "scattering term" does not include any integral over the past history of the system. 
This is equivalent to assume that any collision in the plasma does not depend upon the 
previous ones. On the contrary, quantum distributions posses strong memory effects 
which are relevant for the computation of the final baryon asymmetry since they lead 
to a slowdown of thermalization times and therefore to longer stages of nonequilibrium. 

The paper is organized as follows. In section 2 we give a brief description of the 
basic features of the nonequilibrium quantum field theory and the CTP formalism. 
In sections 3 and 4 we compute the quantum transport equations for local particle 
asymmetries in the bosonic and fermion case, respectively. Sections 5 and 6 contain 
the explicit computation of the CP-violating sources for the right-handed stop and the 
Higgsino numbers. We conclude with an outlook of our findings and comments about 
their implications in section 7. 

2. Some basics of non-equilibrium quantum field theory 

In this section we will briefly present some of the basic features of the nonequilibrium 
quantum field theory. The interested reader is referred to the excellent review by Chou 
et al. J[7) for a more exhaustive discussion. 

The ordinary quantum field theory at finite temperature, which mainly deals with 
transition amplitudes in particle reactions, is not useful to study the dynamics of particle 
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asymmetries. This is because we need their temporal evolution with definite initial 
conditions and not simply the transition amplitude of particle reactions with fixed initial 
and final conditions. 

The most appropriate extension of the field theory to deal with these issues it to 
generalize the time contour of integration to a closed-time path. More precisely, the 
time integration contour is deformed to run from — oo to +00 and back to —00. 

The CTP formalism (often dubbed as in-in formalism) is a powerful Green's func- 
tion formulation for describing nonequilibrium phenomena in field theory. It allows 
to describe phase-transition phenomena and to obtain a self-consistent set of quantum 
Boltzmann equations. The formalism yields various quantum averages of operators eval- 
uated in the in-state without specifying the out-state. On the contrary, the ordinary 
quantum field theory (often dubbed as in-out formalism) yields quantum averages of 
the operators evaluated with an in-state at one end and an out-state at the other. 

The partition function in the in-in formalism for a complex scalar field is defined to 

be 



J, J 1 " 



= Tr T |^exp 
= Tr T + ^exp 
x T_ ( exp 



i J c (j0 + jV f )])p 



+ 44 



(1; 



where the suffic C in the integral denotes that the time integration contour runs from 
minus infinity to plus infinity and then back to minus infinity again. The symbol p 
represents the initial density matrix and the fields are in the Heisenberg picture and 
defined on this closed time contour. 

As with the Euclidean time formulation, scalar (fermionic) fields <fi are still periodic 
(anti-periodic) in time, but with 4>(t,x) = (fi(t — i/3,x), (3 = 1/T. The temperature 
appears due to boundary condition, but now time is explicitly present in the integration 
contour. 

For non-equilibrium phenomena and as a consequence of the time contour, we must 
now identify field variables with arguments on the positive or negative directional 
branches of the time path. This doubling of field variables leads to six different real-time 
propagators on the contour |T7|] . It is possible to employ fewer than six since they are 
not independent, but using six simplifies the notation. For a generic bosonic charged 
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scalar field they are defined as 

G>(x,y) = -i(<f>(x)<f>\y)), 

G<(x,y) = -i(^{y)cj>{x)), 

Gl(x,y) = 9(x,y)G^(x,y)+9(y,x)G^(x,y), 

Gl(x,y) = 9(y,x)G^(x,y)+9{x,y)G^{x,y), 

G r ^(x, y) = G^ — G^ = G^ — G^, G^(x, y) = G^ — G^ = G^ — G^, (2) 

where the last two Green functions are the retarded and advanced Green functions 
respectively and 8(x,y) = 9(t x — t y ) is the step function. For a generic fermion field ijj 
the six different propagators are analogously defined as 

G^(x,y) = -i(ip(x)$(y)), 

G^,(x,y) = +i(ip(y)ilj(x)}, 

G^(x,y) = 9{x,y)G>(x,y) + 9{y,x)G<{x,y), 

G%(x,y) = 9(y,x)G^(x,y) + 9(x,y)G^(x,y), 

Gf( x , y) — G^- G^ = G^ - G^, G^(x, y) = — = G^ - G 1 ^. (3) 

For equilibrium phenomena, the brackets (• • •) imply a thermodynamic average over all 
the possible states of the system. For homogeneous systems in equilibrium, the Green 
functions depend only upon the difference of their arguments (x, y) = {x — y), and 
there is no dependence upon (x + y). For systems out of equilibrium, the definitions 
(0) and (0) have a different meaning. The bracket no longer signifies thermodynamic 
averaging since the concept is now ill-defined. Instead, the bracket means the need to 
average over all the available states of the system for the non-equilibrium distributions. 
Furthermore, the arguments of the Green functions (x, y) are not usually given as the 
difference (x — y). For example, non-equilibrium could be caused by transients which 
make the Green functions depend upon (t x ,t y ) rather than (t x — t y ). 

For interacting systems whether in equilibrium or not, one must define and calculate 
self-energy functions. There are six of them: £*, £*, S < , S > , S r and S a . The same 
relationships exist among them as for the Green functions in (^) and (|3|), such as 

S r = S* - £< = £> - £*, S a = S* - S> = S< - £*. (4) 

The self-energies are incorporated into the Green functions through the use of Dyson's 
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equations. A useful notation may be introduced which expresses four of the six Green 
functions as the elements of two-by-two matrices |L8| 



G ~ [ G> -G l ) ' S ~ { S> -S f J ' (5) 
where the upper signs refer to bosonic case and the lower signs to fermionic case. For 
systems either in equilibrium or non-equilibrium, Dyson's equation is most easily ex- 
pressed by using the matrix notation 

G(x,y) = G°(x,y) + / d 4 x 3 / d 4 x 4 G°(x, x 3 )T,(x 3 , x 4 )G(x 4 , y), (6) 



where the superscript "0" on the Green functions means to use those for noninteracting 
system. This equation is illustrated in Fig. 1, where the thick solid lines represent the 
full Green function and the thin solid lines represent the propagators for the noninteract- 
ing theory. The expression appears quite formidable; however, some simple expressions 
may be obtained for the respective Green functions. It is useful to notice that Dyson's 
equation can be written in an alternate form, instead of @, with G° on the right in the 
interaction terms, see Fig. 2: 

G(x,y) = G°(x,y) + d A x 3 / d A x 4 G(x, x 3 )S(x 3 , x 4 )G°(x 4 , y). (7) 



Eqs. (||) and (0) are the starting points to derive the quantum Boltzmann equations 
describing the temporal evolution of the CP-violating particle density asymmetries. 

3. QBE for particle density asymmetry: the bosonic case 



Kadanoff and Baym |L9[ provided a general method of deriving the QBE's. Here 
we adopt their technique and approach to derive the QBE's for some generic bosonic 
particle asymmetry. This will allow us to derive in a self-consistent way the CP-violating 
sources fueling electroweak baryogenesis in the diffusion equation for the right-handed 
stop asymmetry. 

Our goal is to find the QBE for the following CP-violating current 

( J $( x )) = i(<p\x) d x 4>(x)) = [n<j,(x), Xp(x) 

The zero- component of this current n$ represents the number density of particles minus 
the number density of antiparticles and is therefore the quantity which enter the diffusion 
equations of supersymmetric electroweak baryogenesis. 
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Since the CP-violating current can be expressed in terms of the Green function 



W*)> 



d%-dZ)G<(x,y) 



x=y 



(9) 



the problem is reduced to find the QBE for the Green function (x, y) . To make 
contact with the standard derivation of the QBE , we may go to a center-of-mass 
coordinate system 

X = (T,R) = ±(x + y), (t,r)=x-y. (10) 

Note that T now means the center-of-mass time and not temperature. The notation on 
the Green function is altered to these center-of-mass coordinates 

G<(x, y) = G<(t, f, T, R) = (t -\t,R- V) cj> (t + It, R + if)}. (11) 

The identification x = y in Eq. (|9]) is therefore equivalent to require t — f — 0. 

Our interest is in finding an equation of motion for the interacting Green function 
when the system in not in equilibrium. Such an equation can be found from (|6|) 
by operating by [p x +m 2 ^ on both sides of the equation. Here m represents the mass 
term of the field 0. On the right side, this operator acts only on G°, 

(d x +m 2 ) G^x, y) = S (4 \x, y)h + J d 4 x 3 ^(x, x 3 )G (x 3 , y), (12) 

where I is the identity matrix. It is useful to also have an equation of motion for the 
other variable y. This is obtained from (|7|) by operating by (o y +m 2 ) on both sides of 
the equation. We obtain 



G^x.y) (p y +m 2 ) = 5 (A \x,y)h + J d^G^x.x^H^x^y). 



(13) 



The two equations (|T2|) and (p~3|) are the starting point for the derivation of the QBE for 
the particle asymmetries. Let us extract from ([12]) and (|T3| ) the equations of motions 
for the Green function (x, y) 



+m 2 )G}(x,y) 



d 4 x 3 [^(x,x 3 )G<(x 3 ,y) - Z<(x,x 3 )Gl(x 3 ,y)} , (14) 
d A x 3 \Gl(x,x 3 )E^(x 3 ,y) - G^(x,x 3 )T^(x 3 ,y)\ . (15) 



If we now substract the two equations and make the identification x = y, the left-hand 
side is given by 

drjx) (dn 4 



QX 



dH-d^)G<(x,y) 



x=y 



dX» 



dT 



+ Vr -J 4 



(16) 



and the QBE for the particle density asymmetry is therefore obtained to be 



dn^X) 
dT 



d 4 X3 



G<E* 



(17) 



x=y 



In order to examine the "scattering term" on the right-hand side of Eq. (fLTj) , the first 
step is to restore all the variable arguments. Setting t = f = in the original notation 
of E^(x, x 3 )G ( j } {x 3) y) gives (X, x 3 )(x 3 ,X) for the pair of arguments 
dn^X) 



dT 



d 4 x 3 E* (X, x 3 )G<(x 3 , X) - E<(X, x 3 )Gl (x 3 , X) 



Gl(X, x 3 )E<(x 3 , X) - G<(X, x 3 )E|(x 3 , X) 



The next step is to employ the definitions in (^) to express the time-ordered functions 
G^, E^, and E^ in terms of G5, G^, E^ and G^. Then the time integrals are 
separated into whether t 3 > T or t 3 < T and the right-hand side of Eq. fll8|) reads 



yi>/^< I /^~f<y> yi<^-f> /^f>y< 



= - J d 4 x 3 {0{T-i 
+ Q(t 3 —T) E^G^ + G^E^ — E^G^ — G^E^ |. 

The term with t 3 > T all cancel, leaving T > t 3 . Rearranging these terms gives 



(19) 



dn^X) 
dT 



Vr -MX) 



d 3 r 3 



E>(X, x 3 )G<(x 3 , X) - G>(X, x 3 )E<(x 3 , X) 



+ G<(X, x 3 )X>(x 3 , X) - E<(X, x 3 )G>(x 3 , X) 



(20) 



This equation is the QBE for the particle density asymmetry we were looking for. 
The right-hand side represents the "scattering" term. In the particular case in which 
interactions conserve the number of particles and the latter are neither created nor 
destroyed, their number asymmetry is conserved and should obey the equation of 
continuity dn ( j > /dT+ V_r -j^ = 0. To check that this is indeed the case, one can observe 
that under the assumption that interactions do not change the number of particles, 
most self-energy expressions can be expressed in the following form 



2%(x,y)=g(x,y) G%(x,y) , Ej(ar, y) = g(x, y) G^(x, y) 



(21) 



where g(x, y) = g(y, x) and m is a positive integer. This form of the self-energy is found, 
for instance, for a A |0| 4 theory, where m — 3. In such a case, the terms in the integrand 
of the scattering integral all cancel since 



G>(X,x 3 ) G<(x 3 ,X)- G<(X,x 3 ) G>(x 3 ,X)\x[g(X,x 3 )-g(x 3 ,X)] = 0. 

(22) 



The equation of continuity is therefore satisfied by the QBE. In the most interesting 
cases, however, the particle asymmetries are not conserved in a given environment. This 
occurs if the interactions themselves do not conserve the particle number asymmetries 
and there is some source of CP-violation in the system. Now, if one follows the spirit 
of the usual derivation of Fick's law and the diffusion equation, one should perform 
a simultaneous expansion to first order in the deviations of the distribution function 
n^{X) from its equilibrium value r&, in derivatives of (n^ — nV) and in the particle 
number violating interactions. What is unusual, however, in Eq. ( p0|) is the presence 
of the integral over the time. The physical interpretation of this integral over the 
past history of the system is straightforward: it leads to the typical "memory" effects 
which are observed in quantum transport theory |15], [TJ|. In the classical kinetic theory 
the "scattering term" does not include any integral over the past history of the system 
which is equivalent to assume that any collision in the plasma does not depend upon the 
previous ones. On the contrary, quantum distributions posses strong memory effects and 
the thermalization rate obtained from quantum transport theory may be substantially 
longer than the one obtained from classical kinetic theory. This observation is relevant, 



for instance, when analysing the properties of the quark-gluon plasma pofl . We will 
return to this point in the following. 

The right-hand side of Eq. (^H), through the general form of the self-energy S^, 
contains all the informations necessary to describe the temporal evolution of the par- 
ticle density asymmetries: particle number changing reactions and CP-violating source 
terms, which will pop out from the corresponding self-energy Sep. Notice that so far 
we have not made any approximation and the computation is therefore valid for all 
shapes and sizes of the bubble wall expanding in the thermal bath during a first-order 
electroweak phase transition. If the interactions of the system do not violate CP, there 
will be no CP-violating sources and the final baryon asymmetry produced during su- 
persymmetric baryogenesis will be vanishing. What is noticeable is that we have been 
able to rigously derive a set of quantum transport equations which incorporate the CP- 
violating sources in a self-consistent way. This is an improvement with respect to recent 
treatments where the various CP-violating currents induced by the wall are first derived 
and then converted into sources for the diffusion equations. We will explicitly derive the 
CP-violating source for the right-handed stop number asymmetry and comment about 
its interpretation as a "scattering" term after we have derived the quantum transport 
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equations for fermionic particle number asymmetries. 



4. QBE for particle density asymmetry: the fermionic case 

In this Section we will derive the QBE for the following generic fermionic CP- 
violating current 



(J#(z)) = (ip(x)^^(x)} = \nj,(x),£i,(x) 



(23) 



where ip indicates a Dirac fermion and 7 M represent the usual Dirac matrices. Again, the 
zero-component of this current represents the number density of particles minus the 
number density of antiparticles and is therefore the relevant quantity for the diffusion 
equations of supersymmetric electroweak baryogenesis. 

Our initial goal is to find a couple of equation of motions for the interacting fermionic 
Green function G^,{x,y) when the system is not in equilibrium. Such equations may 
be found by applying the operators (i jd x — M^j and (i f) y +M^j on both sides of Eqs. 
(U) and (0), respectively. Here M represents the mass term of the fermion ip. We find 



i fix -Mj G^(x,y) = 8^{x,y)h + J d 4 x 3 E^{x,x 3 )G^{x 3 ,y), (24) 
Gi,{x,y)(ifl y +M) = -5 {4 \x,y)h- [ d 4 x 3 G4x,x 3 ) Zf(x 3 ,y). (25) 



We can now take the trace over the spinorial indeces of both sides of the equations, sum 
up the two equations above and finally extract the equation of motion for the Green 
function G^ 



Tr 



i fix A 



G %{ x ,v)\ = / d ^sTr E^(x,x 3 )G^(x 3 ,y) -Y>^(x,x 3 )G^,{x 3 ,y) 



- 0^(x,x 3 )Et(x 3 ,y) + G t Jx,x 3 )Y&(x 3 ,y) 



(26) 



If we now make use of the definitions (|Hf), we can work out the left-hand side of Eq. 
(P) 



Tr \i f) x G> (T, R, t,r) + G> (T, R, t, r)i fi 
d 



t=r=0 



d 



t=r=0 



- 7 M 



(27) 
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The next step is to employ the definitions in ([3]) to express the time-ordered functions 
G^p, G^p, E^,, and S^, in terms of G^, G^, and G^. The computation goes along the 
same lines of the analysis made in the previous section and we get 

dn^jX ) - 

or 

+ G<(X, x 3 )£>(x 3 , X) - £<(X, x 3 )G>(x 3 , X) 



+ V R -MX)= I d 3 r 3 f dt 3 Tr \E>(X,x 3 )G<(x 3 ,X) - G>(X,x 3 )E<(x 3 

J J — oo L 



This is the "diffusion" equation describing the temporal evolution of a generic fermionic 
number asymmetry n^. As for the bosonic case, all the informations regarding particle 
number violating interactions and CP- violating sources are stored in the self-energy 
In the following we will explicitly work out the CP-violating sources for charged and 
neutral Higgsinos. 

5. The CP-violating source for the right-handed stop number 

As we mentioned in the introduction, a strongly first order electroweak phase transi- 
tion can be achieved in the presence of a top squark lighter than the top quark [12 



In order to naturally suppress its contribution to the parameter Ap and hence preserve 
a good agreement with the precision measurements at LEP, it should be mainly right 
handed. This can be achieved if the left handed stop soft supersymmetry breaking 
mass fnj is much larger than Mz- Under this assumption, only the right-handed stops 
contribute to the axial stop charge. The right-handed stop current J~ associated to 
the right-handed stop tR is given by 

ja tR = i (t R # t R ) . (29) 

To fix our conventions, let us write the interaction terms among the right-handed stop 
tR, the left-handed stop II and the two neutral Higgses H® 2 , which are responsible for 
the CP-violating source in the diffusion equation for the right-handed stop number nj R 

£ = h t t L {A t H% - fi*H°) t* R + h.c. (30) 

Here the soft trilinear term A t and the supersymmeric mass term fi are meant to be 
complex parameters so that lm(A t fi) is nonvanishing. Even though in this paper we will 
restrict ourselves to the computation of the CP- violating source in the diffusion equation 
of the particle asymmetry nj R , it is clear that the self-energy of the right-handed stop 
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contains the informations about all the other interactions which are responsible for 
changing nj R in the plasma. A typical example is provided by the interaction among 
the right-handed stop, the left-handed top ti and the Higgsino H% . Eq. (p0|) is the QBE 
describing the right-handed stop number asymmetry. Solving this equation represents 
an Herculean task since it is integral and nonlinear. This happens because the self- 
energy functions S > and S < are also functions of the full nonequilibrium Green functions 
of other degrees of freedom of the system. We can make some progress, though. Since 
we know that there is no CP-violating source in the diffusion equation of nq- in absence 
of any Higgs configuration describing the bubble wall profile, we first perform a "Higgs 
insertion expansion" around the symmetric phase (Hf(x)) = Vi(x) = (i = 1,2). 
At the lowest level of perturbation, the interactions ( |30D induce a contribution to the 
self-energy of the form 

Z>cp{x,y) = g C p{x,y)G~' > (x J y), S^ P (x, y) = g C p(x, y)G~' < (x, y), (31) 



where G? > and are now the Green functions for the left-handed stop computed in 



and G£'< 
the unbroken phase and 

9cp(x, y) = h 2 t [A* t v 2 (x) - fwi(x)] [A t v 2 (y) - H*v t (y)] . (32) 

If we now insert the expressions (0) and (|32"D into the diffusion equation (p0|), we get 



dnr 

Qrp ' V H J tR ^t R 



Vr-JT=Si r + ---. (33) 



where 



«%„ = -2i / d 3 r 3 { dt 3 [g C p{X,x 3 ) - g C p{x 3 ,X)] 

J J —oo 



x Im 



-oo 

Gf > (X,x 3 )Gf < (a;3,A) 



+ 



Ah' 2 t I d 3 r 3 [ T dt 3 lm(A t fi) [v 2 (X) Vl (x 3 ) - v 2 (x 3 )v 1 (X)] 



x Im 



G% > (X,x 3 )G^' < (x 3 ,X)\+---. (34) 



where the dots represent the other terms describing the particle number violating in- 
teractions. Sj r is the CP-violating source for the right-handed stop number asym- 
metry. Notice that it vanishes if the relative phase of A t fi is zero and if the ratio 
tan/3 (x) = v 2 (x) /vi(x) is a constant in the interior of the bubble wall. The correspond- 
ing diagram is given in Fig. 3 where the thick dashed line stands for the fact that one has 
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to compute the imaginary part of the diagram. The interpretation of the CP-violating 
source as a "scattering" term is therefore straightforward: the CP-violating source is 
built up when the right-handed stops pass across the wall, they first scatter off the 
wall and are transformed into left-handed stops; the latter subsequently suffer another 
scattering off the wall and are converted again into right-handed stops. If CP-violation 
is taking place in both interactions, a nonvanishing CP-violating source <S~ pops out 
from thermal bath. 

In order to deal with analytic expressions, we can work out the thick wall limit and 
simplify the expressions obtained above by performing a derivative expansion 

*M = E - t^-M x ) W - X "T ■ (as) 

^ n\ d(X») n 

The term with no derivatives vanishes in the expansion (|35|), V2(X)vi(X) — v\{X)v2{X) = 
0, which means that the static term in the derivative expansion does not contribute to 
the source Sj r . For a smooth Higgs profile, the derivatives with respect to the time 
coordinate and n > 1 are associated with higher powers of v^/L^, where v u and 
are the velocity and the width of the bubble wall, respectively. Since the typical time 
scale of the processes giving rise to the source is given by the thermalization time of the 
stops 1/IY, the approximation is good for values of L u T^/v u 3> 1. In other words, this 
expansion is valid only when the mean free path of the stops in the plasma is smaller 
than the scale of variation of the Higgs background determined by the wall thickness, 
L u , and the wall velocity u w . A detailed computation of the thermalization rate of the 
right-handed stop from the imaginary part of the two-point Green function has been 
recently performed in |21[] by making use of improved propagators and including resum- 
mation of hard thermal loopsQ. The thermalization rate has been computed exactly at 
the one-loop level in the high temperature approximation as a function of the plasma 
right-handed stop mass mj R (T) and an estimate for the magnitude of the two- loop con- 
tributions which dominate the rate for small mqr (T) was also given. If mj R (T) > T, 
the thermalization is dictated by the one-loop thermal decay rate which can be larger 
than T [^T[][|. With such value, our derivative expansion is perfectly justified since the 
wall thickness can span the range (10 — 100)/T. 

1 The left-handed stop is usually considered to be much heavier than T and its decay width corre- 
sponds to the one in the present vacuum. 

2 For smaller values of m~^(T), when the thermalization is dominated by two-loop effects (i.e. scat- 
tering), I~ may be as large as 10~ 3 T pi] , 
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The term corresponding to n = 1 in the expansion fl3"5| ) gives a contribution to the 
source proportional to the function 

vi(X)d£v 2 (X) - v 2 (X)d^ Vl (X) = v 2 (X)d%P(X), (36) 

which should vanish smoothly for values of X outside the bubble wall. Here we have 
denoted v 2 = v\ + v\. Since the variation of the Higgs fields is due to the expansion of 
the bubble wall through the thermal bath, the source Sj will be linear in v w . This result 
explicitly shows that we need out of equilibrium conditions to generate the source and 
that we have to call for the CTP formalism to deal with time- dependent phenomena. 
To work out exactly Sj r one should know the exact form of the Green functions which, 
in ultimate analysis, are provided by solving the complete set of Quantum Boltzmann 
equations. However, any departure from thermal equilibrium distribution functions is 
caused at a given point by the passage of the wall and, therefore, is 0{vj). Since 
the source is already linear in v u , working with thermal equilibrium Green functions 
amounts to ignoring terms of higher order in v w . This is as accurate as the bubble wall 
is moving slowly in the plasma. 

The generic finite temperature, real-time propagator in equilibrium G^(k,t x — t y ) 
can be written in terms of the spectral function p<^(k, k ) J22| 

G°/(k, t x -t y ) = / + °° ^e-^-^p^k, k°) { \l + n^k )} 6{t x - t y ) + n^(k°)6{t x - t y )\ 

(37) 

where n^(k°) represents the Bose-Einstein distribution function. 

To account for interactions with the surrounding particles of the thermal bath, par- 
ticles must be substituted by quasiparticles, dressed propagators are to be adopted 
(the use of the full corrected propagators should be done with some care to avoid an 
overcounting of diagrams [p3fl ) and self-energy corrections at one- or two- loops to the 
propagator modify the dispersion relations by introducing a finite width r^(fc). In the 
limit of small decay width, the spectral function is expressed by 

1 1 



(k, k°) 



(3* 



{k° + ie + iT^ - co 2 (k) (k° -%e- iTtf - u%k)_ 
where w^(k) = k 2 + m\(T) and m^(T) is the thermal mass. Performing the integration 
over k° one gets |22| 

G°^(k,t x -t y ) = -J_{[l + n (^-ir )] e -K-W-^) + n(^ + ^)e-^ +ir ^^ 
G$ (k,t x — ty) = (k,t y — t x ). 
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Since the Green functions depend only upon the absolute value of the three-momentum, 



the contribution to the source Sj~ from the n — 1 term in the derivate expansion (|3 
vanishes when we select the space coordinates (// = 1, 2, 3). Indeed, in such a case the 
source is proportional to 



We are therefore left with the expression corresponding to Fig. 3 

S^ r (X) = h\ Im (A/x) v 2 (X)P(X) I 7r , 

where $(X) = df3(X)/dt x , 

k 2 



(40) 



(41) 



X: 



dk 



2n 2 UJr OJr 



X 



1 + 2Re(n r j) 7(u^, F 7r , ^, r~ ) + (l + 2Re(n ? j) I^F 
2 (im(n^) + Im(n r j) G(u^ R , F Tr ,u^ l ,F T[ 



(42) 



and n 



l/[exp (^ {L) /T + zr ?ii(L) 
J(a, 6, c, 6?) : 



/TJ — 1 . The functions / and G are given by 

a + c 



+ 



G(a, b, c, g!) 



1 


1 


2 [(a + c) 


2 + (6 + rf) 2 ] 


1 


1 


2 [(a - c) 


2 + (6 + d) 2 ] 


1 


1 


2 [(a + c) 


2 + (6 + rf) 2 ] 


1 


1 



sin 



sm 



cos 



COS 



2arctan 
2arctan 
2arctan 
2arctan 



b + d. 
a — c 
b + d. 
a + c 



6 + 
a — c 



6 + 6? 



(43) 



co>r . This resonant 



2 [(a - c) 2 + (6 + df] 

Notice that the function G(uj r , Fj r , Fj) has a peak for ujj r 
behaviour || is associated to the fact that the Higgs background is carrying a very low 
momentum (of order of the inverse of the bubble wall width L^) and to the possibility 
of absorption or emission of Higgs quanta by the propagating supersymmetric particles. 
The resonance can only take place when the left-handed stop and the right-handed 
stop do not differ too much in mass. By using the Uncertainty Principle, it is easy to 
understand that the width of this resonance is expected to be proportional to the ther- 
malization rate of the particles giving rise to the baryon asymmetry. Within the MSSM, 
however, it is assumed that m 7 ^> T and the resonance can only happen for momenta 
larger than mj . Such configurations are exponentially suppressed and do not give any 
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relevant contribution to the CP- violating source. Nonertheless, if the electroweak phase 
transition is enhanced by the presence of some new degrees of freedom beyond the ones 
contained in the MSSM, e.g. some extra standard model gauge singlets, the resonance 
effects in the stop sector might be relevant. What is relevant here is that the source 
may be dominated by particles with long wavelengths in direction perpendicular to the 
wall for which the classical approximation breaks down. 

6. The CP-violating source for the Higgsino number 

The Higgs fermion current associated with neutral and charged Higgsinos can be 
written as 

Jt = HYH (44) 

where H is the Dirac spinor 

ff-(f) m 

and H 2 = H 2 (H}), Hi = H® (iff) for neutral (charged) Higgsinos. The interac- 
tions among the charginos and the charged Higgsinos which are responsible for the 
CP-violating source in the diffusion equation for the Higgs fermion number read 



£ = -g 2 {if [v x {x)P L + e id »v 2 (x)P R ] w} + h.c, (46) 

where 9^ is the phase of the /i-parameter and Pl,r are the chirality projector operators. 
Analogously the interactions among the Bino, the W^-ino and the neutral Higgsinos 
are 

C = -^H^[v 1 (x)P L + e w ^v 2 (x)P R ] (g 2 W 3 -9iB)}+ h.c. (47) 

To compute the source for the Higgs fermion number we again perform a "Higgs 
insertion expansion" around the symmetric phase. At the lowest level of perturbation, 
the interactions of the charged Higgsino induce a contribution to the self-energy of the 
form (and analogously for the other component S^ P ) 



S<p(x, y) = gZ P (x, y)P L G°£(x, y)P L + ggp(x, y)P B G°^(x, y)P R , (48) 



where 



9cp( x ,v) = 9 2 2 vi{x)v 2 {y)e *», 

9cp(x,y) = g 2 My)v 2 (x)e ie ». (49) 
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Similar formulae hold for the neutral Higgsinos. 

Analogously to the case of right-handed stops, the dispersion relations of charginos 
and neutralinos are changed by high temperature corrections [^J . Even though fermionic 
dispersion relations are highly nontrivial, especially when dealing with Majorana fermions 
relatively simple expressions for the equilibrium fermionic spectral functions may 



be given in the limit in which the damping rate is smaller than the typical self-energy 
of the fermionic excitation |TB[. For instance, the spectral function of the charged 
Higgsinos may be written as 



(50) 



(fc° + is + %V~) 2 - u 2 ~{k) (k° -is- iT~) 2 - u 2 ~(k) 



where uJ 2 ~(k) = k 2 + m~(T) and m~(T) is the Higgsino effective plasma squared mass in 
the thermal bath which may be well approximated by its value in the present vacuum, 
m~(T) ~ | /x | 2 . Similarly, should be replaced by M 2 for p^(k,k°), and by M 1 for 
pg(k, A; ). Inserting the expressions fl48|) and (|49|) into the diffusion equation (p8|) , we 
can perform a Higgs insertion expansion of the CP-violating source. The computation 
goes along the same lines of the calculation done in the previous section and it is easy 
to show that the CP- violating source 



H 



d d r 3 



dU Tr 



— oo 

> 



E> P (X, x 3 )G|(x 3 , X) - G|(X, x 3 )S<p(x 3 , X) 



+ Gf (X, x 3 )S> P (x 3 , X) - £< P (X, x 3 )G|(x 3 , X) 
containes in the integrand the following function 



(51) 



c/ CP (X,x 3 )+^cp(^,a;3)-^cp( a; 3,^)-^cp( a:; 3,^) = 2isin^ [v 2 {X)v 1 {x z ) - v 1 (X)v 2 (x 3 )} 

(52) 

which vanishes if Im(/i) = and if the tan (5{x) is a constant along the Higgs profile. 
Performing the "Higgs derivative expansion", we finally get 



S S {X) = Inn» v\X)(3{X) 



3M 2 g 2 2 l^ + M l g 2 1 l ] 



H 



(53) 



where 



1 



w 

H 



oo h% 

dk— 
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(l - 2Re(n^)) I(w s , V s , u^, T~) + (l - 2Re(n 5 )) T^, r~) 

+ 2(lm(n g ) + lm(n^))G(u s ,T s ,u^,T^) 



and w l(vi>) = k 2 + \v\ 2 (Mi) while n g(#) = 1/ [exp (u S( ^/T + iT S( ^/T) + 1 



'H(W)l J - 1 l,J -fl-(W^)/ 



(54) 



The 



exact computation of the damping rate of charged and neutral Higgsinos will be given 
elsewhere [|2lj|| . The Bino contribution may be obtained from the above expressions by 
replacing M2 by M\. As for , the CP- violating source for the Higgs fermion number is 
enhanced if M2, Mi ~ [i and low momentum particles are transmitted over the distance 
L w . This means that the classical approximation is not entirely adequate to describe the 
quantum interference nature of CP- violation and only a quantum approach is suitable 
for the computation of the building up of the CP-violating sources. 

7. Outlook Let us now look back and comment about the various aspects of our 
findings. 

-Comparison to previous work- 
One of the merits of the CTP formalism is to guide us towards a rigorous and 
self-consistent definition of the CP-violating sources within the quantum Boltzmann 
equations. On the contrary, previous treatements |S| are characterized by the fol- 
lowing common feature: CP-violating currents were first derived and then coverted 
into sources for the diffusion equations. This procedure is (at least theoretically) not 
self-consistent. More specifically, CP-violating sources S associated to a generic charge 
density j° were constructed from the current by the definition S = doj° [0, ||. A rig- 
orous computation of the CP-violating currents for the right-handed stop and higgsino 
local densities was performed in || by means of the CTP formalism. Since currents 
were proportional to (3 in the tick bubble wall limit, i.e. proportional to the first time 
derivative of the the Higgs profile, sources turned out to be proportional to the second 
time derivative of the Higgs profile ||. Our results, however, indicate that the sources 
in the quantum diffusion equations are proportional to the first time derivative of the 
Higgs configuration. A comparison between the sources Sj r , see Eq. (0), and S^, see 
Eq. (|53D, obtained in the present work and the currents j given in Eqs. (14) and (21) 
of ref. 0] indicate that they may be related as 

S(T) ~ (55) 
r 
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and it may be interpreted as the time derivative of the current density accumulated at 
the time T after the wall has deposited at a given specific point the current density j° 
each interval r 

S(T) ~d J T dt 3 -^-. (56) 

Here r = T" 1 is the thermalization time of the right-handed stops and higgsinos, re- 
spectively. The integral over time is peculiar of the quantum approach and it induces 
memory effects. This tells us that the source obtained self-consistently in the present 
work differs from the one adopted in [BJ by a factor ~ L^T/v^ (in the rest frame of the 
advancing bubble wall). Since L^T/v^ > 1 for the Higgs derivative expansion to hold, 
this result is important as far as the numerical estimate of the final baryon number is 
concerned. We notice that the definition of the source given in [[7| is very similar to 
(p5|) even though it was not motivated by first principles and it did not incorporate 
self-consistently the decoherence effects which have a crucial impact on the generation 
of the CP-violating observables. 
-Memory effects- 

Baryogenesis is fueled when transport properties allow the CP-violating charges to 
efficiently diffuse in front of the advancing bubble wall where anomalous electroweak 
baryon violating processes are not suppressed. However, the CP-violating processes of 
quantum interference, which build up CP-violating sources, must act in opposition to 
the incoherent nature of plasma physics responsible for the loss of quantistic interference. 
If the particles involved in the process of baryon number generation thermalize rapidly, 
CP-violating sources loose their coherence and are diminished. The CTP formalism 
properly describes the quantum nature of CP-violation and tells us that CP-violating 
sources evaluated at some time T are always proportional to an integral over the past 
history of the system. Therefore, it is fair to argue that these memory effects lead 
to "relaxation" times for the CP-violating sources which are typically longer than the 
ones dictated by the thermalization rates of the particles in the thermal bath. In 
fact, this observation is valid for all the processes described by the "scattering" term 
in the right-handed side of the quantum diffusion equations. The slowdown of the 
relaxation processes may help to keep the system out of equilibrium for longer times 
and therefore enhance the final baryon asymmetry. There are two more reasons why 
one should expect quantum relaxation times to be longer than the ones predicted by 
the classical approach. First, the decay of the Green's functions as functions of the 
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difference of the time arguments: an exponential decay is found in thermal equilibrium 
when one ignore the frequency dependence of self-energies in the spectral functions, 
e.g. |G > (k,t,t , )| ~ l^k)] x exp [-r(k,u)\t - t'\}. The decay of the Green's functions 
restrict the range of the time integration for the scattering term, reduces the integrals 
and, therefore, the change of the local particle number densities as a function of time. 
The second effect is the rather different oscillatory behaviour of the functions G > and 
G < for a given momentum, as functions of the time argument difference. 
-Resonance effects- 

In the limit of thick bubble walls, the CP-violating sources are characterized by 
resonance effects [|J when the particles involved in the construction of the source are 
degenerate in mass. The resonance is manifest in the function G defined in fl43l). The 
interpretation of the resonance is rather straightforward if we think in terms of scatter- 
ings of the quasiparticles off the advancing low momentum bubble wall configuration. 



A similar effect has been found in ref. [|K| where the system was studied in the clas- 
sical limit. These classical treatments somehow obscure the origin of the CP- violating 
effects as resulting from quantum interference and the origin of the resonance is less 
transparent to us. However, these methods should provide reasonable approximations 
to our formulae to those particles whose wavelength is short compared to v^/T . On the 
other hand, formulae should not agree for small V because our source is dominated by 
particles with long wavelength. In this regime, the classical approximation breaks down 
since it requires that the mean free path should be larger than the Compton wavelength 
of the underlying particle. This is relevant because quasiparticles with long wavelengths 
give a significant contribution to CP-violating sources. 
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Figure Captions 



Fig. 1 : The first Dyson equation for the Green function matrix G. Here thick and 
thin solid lines represent the Green functions for the fully interacting system and for 
the free theory, respectively. 

Fig. 2 : The second Dyson equation for the Green function matrix G. The meaning 
of thick and thin solid lines is as in Fig. 2. 

Fig. 3 : The Feynman diagram representing the CP-violating source for the right- 
handed stop number. The indeces i,j run from 1 to 2 and the solid dashed line means 
that the imaginary part of the diagram should be considered. 

Fig. 4 : The Feynman diagram representing the CP-violating source for the Hig- 
gsino number. The indeces i,j run from 1 to 2 and the solid dashed line means that 
the imaginary part of the diagram should be considered. 
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